CALIBRATION OF CHAOTIC MODELS FOR INTEREST RATES 
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Abstract. In this paper we calibrate chaotic models for interest rates to market data using a 
polynomial-exponential parametrization for the chaos coefficients. We identify a subclass of one- 
variable models that allow us to introduce complexity from higher order chaos in a controlled way 
while retaining considerable analytic tractability. In particular we derive explicit expressions for 
bond and option prices in a one-variable third chaos model in terms of elementary combinations of 
normal density and cumulative distribution functions. We then compare the calibration performance 
of chaos models with that of well— known benchmark models. For term structure calibration we find 
that chaos models are comparable to the Svensson model, with the advantage of guaranteed positivity 
and consistency with a dynamic stochastic evolution of interest rates. For calibration to option data, 
chaos models outperform the Hull and White and rational lognormal models and are comparable to 
LIBOR market models. 
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1. Introduction. The purpose of this paper is to investigate the calibration 
performance of interest rate models based on the Wiener chaos expansion. Chaotic 
models were introduced in [16| as an axiomatic framework for interest rates satisfying 
both no arbitrage and positivity conditions, following the line of research initiated by 

, where zero-coupon bond prices are modeled as 

M ts ds 

for a parametrized family of positive martingales M ts . It was then shown in [25] and 
|26j that one can focus instead on modeling a supermartingale V t which is related to 
the martingales M ts through 

/oo 
M ts ds, (1.2) 

The chaotic approach is derived from the observation that Vt can itself be written as 
the conditional variance of a terminal random variable X^. This square integrable 
random variable then has a unique Wiener chaos decomposition, and studying the 
way its expansion coefficients affect the corresponding interest rates model becomes 
the subject of the theory. 

Notice the successive simplifications in the class of objects that one needs to 
model: from an entire family of martingales M ts in the Flesaker-Hughston framework, 
to the process Vt in the potential approach, and finally the single random variable 

in chaotic models. Arguably, each step of the way reduces the arbitrariness in the 
modeling exercise, with presumed advantages for calibration to real data. In particu- 
lar, the chaos expansion allows one to successively introduce randomness in the model 
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in a way that can, in principle, be made to correspond to the increased complexity of 
financial instruments under consideration. Such are the types of statements that we 
propose to put to test in this paper. 

Models in the Flesaker-Hughston framework have been implemented and cali- 
brated, for example, in [5] and [14], whereas implementations of the potential ap- 
proach can be found in |19j and |25j . To our knowledge, we present here the first 
practical implementation and calibration to market data of chaotic interest rate mod- 
els. We adopt the day-by-day calibration methodology used in [TJ], [TS], and [2"5] . 
The motivation for this is to capture the prices of liquid interest rate derivatives such 
as caps and swaptions by a model as parsimonious as possible, which can then be 
used for pricing and hedging of exotic options, such as the Chooser flexible cap and 
Bermudan Swaption. 

After reviewing the chaotic approach in Section [2] we implement to two separate 
calibration exercises. In Section [3] we consider the calibration of chaotic models to 
the observed term structure of interest rates. For comparison, we use the Nelson- 
Siegel [2TJ and Svensson [28] models for forward rates as benchmarks. These so 
called descriptive models [4] are examples of a general exponential-polynomial class 
of models analyzed in pQ, which we use as motivation for the parametric form we 
adopt for the chaos coefficient functions. We find in this section that chaotic models 
perform comparably to descriptive models with the same number of parameters, with 
the advantage of avoiding problems with positivity and consistency. 

We then move to a full calibration to yields and option prices in Section |4j We 
recall known expressions for option prices in a second chaos model, derive the corre- 
sponding formulas in a third chaos model, and calibrate them to two separate data 
sets in three different ways. For comparison, we also calibrate the Hull and White 
model [T7], the rational lognormal proposed in [TT] and a lognormal LIBOR market 
model [12] . We find that chaos models generally perform much better than the Hull 
and White and rational lognormal models and have fitting errors comparable to those 
obtained with a LIBOR model. When we apply an information criterion that takes 
into account the different number of parameters in the models, we find that one of 
our third chaos models consistently outperforms the LIBOR market model. Section 
[5] concludes the paper by summarizing the results and pointing to future research in 



2. The chaotic approach. We review the framework proposed in [IB], whereby 
a general positive interest rate model with no arbitrage opportunities is associated 
with a square integrable terminal random variable X x , which can then be modeled 
using a Wiener chaos expansion. 

2.1. Basic definitions. Let (fi, J 7 , P) be a probability space equipped with the 
standard augmented filtration (J-t)o<t<oo generated by a fc-dimensional Brownian 
motion Wt and suppose that £ L 2 (ft, P) is a random variable with the property 
that is not T t -measurable for any finite value of t and E^X^] = 0. Using the 
martingale representation theorem, we can write 



Jo 

for an adapted process o t . Denoting the conditional variance of X^ with respect to 
the cr-algebra F t by 



the area. 




(2.1) 



V t = E[{X 00 -E t [X 00 ]f\F t ] , 



(2.2) 
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it follows from the Ito isometry that 
V t = E 

Defining the increasing adapted process 

A t = 

we have that 









' POO 








= E 


J °lds 


Ft 













l ds. 



V t = E[A 00 -A t \F t ], 



(2.3) 



(2.4) 



(2.5) 



from which it follows that V% is a potential, that is, a positive supermartingale satis- 
fying E[Vt] — > as t — > oo. It is then well-known (see, for example, [23]) that V t can 
be used as a state-price density for an arbitrage-free interest rate model in which the 
price of a zero-coupon bond with maturity T is given by 



Pit = 



E{V T \Ft} _ Z tT 



Vi 



where 



ZtT — E{Vr\J-t_ 



< t < T < co, 



< t < T < co. 



(2.6) 



(2.7) 



It follows from (2.51 that for each fixed t the processes ZtT, and consequently the 
bond prices PtT, are decreasing functions of the maturity T, which in turn implies 
that all instantaneous forward rates 



ftT = -J^log/Vr 



(2.8) 



are automatically positive in this framework. Furthermore, the short-rate r t and the 
market price of risk vector A* are determined by the dynamics of Vt as follows: 



dV t = -r t Vtdt - XtVtdWt, V > 0. 
For later use, we introduced the family of positive martingales 



M ts = E[al\Ft], 



< t < s < co, 



(2.9) 



(2.10) 



so that Vt can be written as in ( 1.2 1 and bond prices as in ( 1.1 ). 

Returning to the squared integrable random variable X^, the Wiener chaos de- 
composition (see [H] and [3T]) says that it can be written as the L 2 -convergent sum 



Ms)dW s 



o Jo 



cj> 2 {s lSl )dW Sl dW s 



4>z(s, sx,s 2 )dW S2 dW ai dW s + ■ 



(2.11) 



for deterministic square-integrable functions 0i(s), 4>2(s, si), <p3{s, Si, S2), . . . called 
the chaos coefficients of X^. We say that an interest rate model is an n-th chaos 
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model if the decomposition for the random variable can be completely determined 
by its first n coefficient functions. 

The basic insight of the chaotic approach is that the decomposition above pro- 
vides a way to add complexity to an interest rate model in a controlled manner. For 
example, as we will see in the next section, the first order chaos correspond to deter- 
ministic interest rate models, whereas the second order chaos give rise to stochastic 
interest rate models with randomness governed by a parametric family of Gaussian 
processes. 

More importantly for calibration purposes, the increased complexity in the inter- 
est models should be related to the instruments that are available in the market. In 
what follows, we propose a systematic way to calibrate chaotic interest models start- 
ing with bond prices alone and then gradually increasing the complexity of market 
instruments included in the calibration. Our general strategy will consist of choosing 
general parametric forms for the deterministic functions fa, i = 1,2,... and then fit 
the parameters to market data according to the bond and option pricing formulas 
emerging from the models. 

2.2. Option Pricing. We now describe how the prices of the most common 
interest rate options can be written in terms of the quantities defined in Section |2.1| 
Since Vt is a state-price density, it follows that the price at time t of a derivative with 
payoff Ht is given by 

E[V T H T \F t ] 



H t = 



Vi 



(2.12) 



We see that expression (2.6 1 for bond prices is simply a special case of (2.12) with 
Ht = Ptt = 1- We apply this general expression to European put options, caplets 
and swaptions, since these will be the only ones needed in our calibration section, 
although similarly results also hold for other interest rate derivatives, notably call 
options and floors (see [30]). We follow the standard definitions and notations in [2]. 

A put option with maturity t and strike price K written on a bond with maturity 
T >t correspond to the following payoff 



(k - p tT y 



(2.13) 



Using (2.12), we see that the price of the a bond option at time s < t < T is given by 

E[ V t {K - P tT )+\F s ] = E[(KZ tt - Z tT y 

V. v s 



p(s,t,T,K) = 



\?s] 



(2.14) 



A caplet is a call option with strike K and maturity T written on a spot LIBOR 
rate L(t,T) for the time interval [t, T], that is, it is defined by the payoff 

Af(T-t)(L{t,T)-K)+, (2.15) 

where Af is a fixed notional amount. Recalling that the spot LIBOR rate is defined 
as 



we see from ( 2.12[ ) that the value of a caplet at time s < t is 

E [V T N{T - t){L(t,T) - K)+\F S ] 



Cpl(s,t,T,Af,K) 



V, 



E 



(Z tt -Z tT (l + K(T-t)))' 
4 



(2.16) 



(2.17) 



Comparing (2.17) and (2.141 it is easy to deduce the following well-known relation 
H page 41]: 



Cpl(s, t, T,J\[, K) = Af(l + K(T - t))p s, t, T, 



l + K{T-t) 



(2.18) 



which allows us to obtain prices for caplets provided we know how to calculate the 
prices of bond put options. 

Observe that the payoff (2.15) can be rewritten as 

Af(T-t)(F(t,t,T)-K)+, 

where the forward LIBOR rate F(s, t, T) is defined as 



(2.19) 



F(s,t,T) 



1 



P,, 



1 



0<s<i<T<oo. 



T-t \PsT 

Using this, we can define the caplet implied volatility as the value Vt.T satisfying 



(2.20) 



Cpl(0, t, T,Af, K) = JVP 0T {T - i)Black f K, F(0, t, T),v t , T Vi 



where 



Black(if, F, v) := F$ 



\og(F/K) 



log(F/K) 



(2.21) 



(2.22) 



is the Black formula and $(•) denotes the standard normal cumulative distribution 
function. We say that a caplet is at-the-money (ATM) at time t — if K — F(Q, t, T). 
Given a set of maturities Ti, the so-called term structure of caplet volatility is the 
implied volatility curve obtained from ATM caplets of consecutive maturities, that is, 
a function of the form — > vx i _ lt T i - 

A payer Swaption with maturity t is defined by the payoff 



A^(T ( - Ti-i)PtTi {F(t, Tj_i, T,-) - K) 



(2.23) 



where T = (To, T\, . . . , T n ) is a set of future dates with Tq = t and AT is the notional. 
The total length of the time interval (T n — To) is called the tenor of the swaption. 
Recalling the definition of the forward LIBOR rate in (2.20), we obtain from 



( 2.12|) that the value of the payer swaption at time s < t is given by 

V t \Mj2( Ti - T l - 1 )P tTl (F(t,T l _ 1 ,T i ) - K) 

\ i=l 

lz tt -Z tTn -Kf2( T i-Ti-i)Z m ) .F, 



Swp(s,T,AA,A') = W E 



4b 



i=l 



F L 



(2.24) 



Observe that the payoff in (2.23) can be rewritten as 



A/'5^(2;-T < _i)p m ( 1 s(t,i,T n )-jr)J , 
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(2.25) 



i=l 



where the forward swap rate S(s,t,T n ) is defined as 

S(s,t,T n ):= ** S '~% T " , 0< S <<<T„<oo. (2.26) 

Similarly to caplets, we can use this to define the swaption implied volatility as the 
value Ut,T„ satisfying 

n 

Swp(0, r,7V, if) =^/'53(T i - T^^PoTiBlack (if, 5(0, t, T n ),«t, rn Vt) • (2.27) 

i=l 

The swaption is said to be at-the-money (ATM) at time t = if K = S(0,t,T n ). 

This concludes our general overview of chaotic models. We will have a chance to 
explore these formulas in more detail after we introduce specific models of different 
chaos orders in the calibration sections that come next. For simplicity, we will assume 
from now on that the underlying Brownian motion W t is one dimensional, although 
extensions to the multidimensional case should be relatively straightforward. 

3. Term structure calibration. We first consider the calibration of the chaos 
models to the bond prices Pqt observed at time t = 0, or equivalently, to the corre- 
sponding forward rates 

for = -^ lo S p oT- 

Historically, the problem of characterizing the term structure of interest rates 
implied by bond prices at time t — has been extensively tackle by so-called descrip- 
tive models, which consists of choosing a parametric form for the forward rates /qt 
and then fit the relatively small number of parameters to the observed yield curve. 
Popular examples of this approach are: 

( b + [bi + b 2 T]e- ClT (Nelson and Siegcl) 
/ot=< b + [b x + b 2 T\e-^ T + b 3 Te- c ? T (Svensson) (3.1) 
I b o + J2t=i bie~ CzT (Cairns), 

introduced respectively in jjl], |5S] and All of these examples fall within the 
exponential-polynomial family of forward rates analyzed in [1]: 

n hi 

foT = L Q (T)+Y / HT)e- c ' T where L ( (T) = £ fcyT*. (3.2) 

i=l j=0 

As explained in [4], descriptive models provide a snapshot, but not a movie. 
That is to say, they do not provide a model for the stochastic evolution of interest 
rates. Nevertheless, they can be used to price simple instruments such as forward rate 
agreements and swaps, which do not depend on the dynamics of interest rates. They 
have also been widely used by central banks for the purposes of monetary policy (see 
the table presented in [TO]). 

For more complicated interest rate derivatives, it becomes necessary to embed 
such static descriptions into a fully dynamic model. For example, one can use the 
fitted forward rate curve /ot as the initial term structure for a model within the 
Heath- Jarrow-Morton framework |15j . supplementing it with the specification of the 
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volatility structure. In other words, the specification of the initial term structure used 
as an input is done separately from its stochastic evolution. 

By contrast, instead of focusing on the forward rates fox, we propose to use the 
observed term structure at time t = to calibrate the deterministic functions (j>i{T), 
which then completely specify the stochastic evolution of interest rates. Since these 
are markedly different procedures, their performance need to be compared using real 
observed data, as we do later in this section. Before that we need to describe the term 
structure of interest rates for chaotic models of different orders in mode detail. 

3.1. Term structure in first chaos models. In a first chaos model we have 

Xoo = / Ms)dW s (3.3) 



where 4>i(s) is a deterministic function of one variable. Comparing this with (2.1 1 we 
conclude that a s = <fii(s), so that 



V t = E 

which implies that 



<7„ds 



/OO 
Mrfds, (3.4) 



dV t = -<j> x {t) 2 dt. 



Comparing this with (2.9| then gives 



n ~ v t -jrmds (3 - 5) 

X t = 0. (3.6) 

In other words, in first chaos models the state-price density is a deterministic function, 
implying that the interest rates associated with it are themselves deterministic. In 
particular, M ts — 4>\{s) constitutes a family of constant martingales for each value of 
s and bond prices are given by 

E[V T \F t ] flMtsds £^{s}ds _ fTrada 

Vt f t °°M ts ds f t °° <t>l(s)ds ■ [6J) 

It follows that instantaneous forward rates in first chaos models have the form 

which are manifestly positive, as they should be in any chaotic model. Moreover, 
because they represent market expectations at time t for what the (deterministic) 
short-rate will be at time T, instantaneous forward rates in first chaos models are 
independent of t and always equal to tt- 

As already noticed in [TB], these observations show that first chaos models char- 
acterize general arbitrage-free, positive, deterministic interest rate term structures, 
which should not be immediately disregarded as trivial, since they constitute the 
starting point for the majority of the applications of interest rate theory. Chaotic 
models then provide a well-defined way to embed such deterministic models into their 
stochastic generalizations, as we will demonstrate in the remainder of this paper. 
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From the point of view of calibration, we can calibrate a first chaos model to the 
term structure at time t = by focusing on the equation 



Pi 



OT 



J™ (j>i(s)ds 



(3.9) 



f£°M a ) ds ' 

Because of the deterministic nature of first chaos models, this is enough to characterize 
bond prices at all times, since it follows from (3.7) that P t r = Pot /Pot- 



3.2. Term structure in second chaos models. In a second chaos model we 



have 



Ms)dW s 



o Jo 



<p 2 (s,u)dW u dW s 



for deterministic functions </>i(s) and 4>2{s,u). Comparing this with (2.1) gives 



er t = 0i(i)+ / (t> 2 (s,u)dW u . 
Jo 

It then follows from the conditional Ito isometry that 
M ts = E[a 2 s \T t ] 

<px(s)+ cj> 2 {s,u)dW u ) + (j> 2 2 (s,u)du 



(3.10) 



— ^ts Qts ~i~ Qss 

where, for each fixed value s £ [t, oo), the process 

R ts ■= <M S ) + / 02(s,u)dW u 
Jo 

is a martingale in the variable t with quadratic variation 



ts := / 4> 2 {s,u)du. 



(3.11) 



(3.12) 



(3.13) 



Setting t — in (3.11) and using (1.1) we have that the initial term structure of 



bond prices in second chaos models is generically given by 

It (^iO) + It u)du) ds 



Pot 



ir (4>i( s ) + lo 4>l{s, u)du) ds ' 



(3.14) 



Expression (3.14) looks uncomfortably complicated for the purpose of calibrating 
a relatively simple object such as Pot- To achieve greater tractability, we now recall 
the class of factorizable second chaos models introduced in [TB] , whereby 



<^i(s) = a(s), 4>2{s,u) = f3(s)j(u). 



(3.15) 



for deterministic functions a, f3 and 7. In other words, factorizable second chaos 
models correspond to random variables Xoc of the form 



Xcx, — 



a{s)dW s 



Jo 



p(s)j(u)dW u dW s 



(3.16) 



for < u < s < oo. In this case we have that 
R ts = a(s)+P(s)Rt, 

where 



}ts=P 2 {s)Q u 



Rt 



l{u)dW u 



(3.17) 



(3.18) 



is a Gaussian martingale with quadratic variation Q t := J Q ^y 2 (u)du. Inserting (3.17) 
into (3.11) we obtain 



M ts ds 

(a(s) + P{s)R t f - fiQ t + /3 2 (s)Q s 
A(T) + B(T)R t + C(T)(R 2 t - Q t ) 



ds 



(3.19) 



where 



A(T) 



[a 2 (s)+/3 2 (s)Q s }d s 



B(T) = / 2a(s)f3(s)ds 
Jt 



(3.20) 



C{T) 



f3 2 (s)ds. 



Using the fact that Vt — Ztt> it is easy to see that the state-price density in this 
case is given by 



V t = A(t) + B{t)R t + C(t)(R 2 - Q t ), 



(3.21) 



which is a simple quadratic polynomial in the Gaussian process Rt- Therefore bond 
prices in a factorizable second chaos model are given by the following ration of 
quadratic polynomials in R t : 



tr 



A(T) + B(T)R t +C(T)(R 2 -Q t ) 
A(t) + B(t)R t + C(t){R 2 - Q t ) 



In particular, the initial term structure is given by 

A(T) = J^[a 2 ( S ) + P 2 ( S )Q s }d S 
A(0) ~ J^[a 2 ( S ) + P 2 (s)Q s ]d s - 



Pot = 



(3.22) 



(3.23) 



Despite being much simpler than ( |3.14 ), we observe that (3.23) still offers too 
much freedom for the choice of parametric forms for the deterministic functions in- 
volved. To make this point more explicitly, observe that from the strict point of view 
of calibrating the initial term structure, we see that a factorizable second chaos model 
behaves exactly like a hrst chaos model if we identify 



Ms) = a 2 (s)+(3 2 ( S ) / 7 2 ( U )d U) 
Jo 
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(3.24) 



for arbitrary choices of a, j3 and 7. This is to be expected in all chaos models, since the 
initial bond prices available for calibration correspond to the one-variable function 
Pot, regardless of how many deterministic functions we wish to calibrate. Of 
course, once calibrated, chaos models of different orders give rise to very different 
future term structures: a deterministic one with P t T = Pot /Pot for first chaos models 



and a fully stochastic one where PtT is given by (3.22) for factorizable second chaos 
models, for example. 

These remarks point to a general yet delicate feature of chaotic models: one can 
achieve increased stochastic complexity in the form of higher order chaos at the ex- 
pense of having to calibrate complicated expressions that are not necessarily supported 
by the data available. Our approach to deal with this problem consists of increasing 
the complexity of the models as gradually as possible, while carefully evaluating their 
statistical performance when calibrated to real data. 

Accordingly, we propose to investigate the subclass of factorizable second chaos 
models with 7(11) = 1 and /3(s) satisfying 



(3.25) 



f3(s) 2 sds < 00, 

that is, models corresponding to the terminal random variable 

/>00 />OG pS />00 />00 

X x = / a{s)dW s + / (3(s) / dW u dW s = / a{s)dW s + / p{s)W s dW s . (3.26) 
Jo Jo Jo Jo Jo 

As we can see, these one-variable second chaos models provide the simplest stochastic 
extensions of first order chaos, for which we have that Rt = Wt and the initial term 
structure of bond prices is given by 



Pqt — 



J T °[a 2 ( S ) + f3 2 {s)s)ds 
J °°[a^(s)+j3^(s)s]ds' 



(3.27) 



3.3. Term structure in higher order chaos models. In a general chaos 
model wc that 



<M«i) + / 4>2(si,s 2 )dW S2 




/• 00 

Zot= / E 
Jt 

4>i(si) + J 4>l{si,s 2 )ds 2 + 



dsi 



dsi. 



Therefore 



where 



P 



0T 



Zqt = It Hs)ds 
Zot / t °° ^p(s)ds' 



(3.28) 



^i( s ) + /o 4>l{s,u)du 



(first chaos) 
(second chaos) 



^ (S) ~ I 4>\{s) + fa </>|(s, u)du + J Q S f " 0§(s, u, v)dvdu (third chaos) (3 ' 29) 
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In other words, as mentioned in the previous section, equation (3.28) shows that 
from the strict point of view of calibrating the initial term structure Pot, ah chaos 
models behave like first order chaos models, albeit with a complicated choice of func- 
tion ij}(s). The differences between chaos models will become apparent, however, when 
we tackle the calibration to option prices in Section |4j For now, we find unnecessary 
to consider fully general higher order chaos and restrict ourselves to the simplest pos- 
sible third chaos models. That is to say, we start with factorizable third chaos models 
of the form 



a{s)dW s 



OO p S pU 



P(sh(u)dW u dW s 



6(s)e(u)((v)dW v dW u dW s , 



for < v < u < s < oo. We then take 7 = e = £ = 1, j3 satisfying (3.251 and 5 
satisfying 



/•OO 

/ 6(s) 2 s 2 ds < 00 
Jo 



to obtain the class of one-variable third chaos models of the form 



— 



a(s)dW s 



OO ps pU 



P{s)dW u dW s 



Jo 



5{s)dW v dW u dW s 



a(s) + p{s)W s + -5{s){W 2 s 



J0J0 
dW s . 



We see from (3.31) that in this case 

a s = a(s) + (5{s)W s + ^S(s)(W 2 - s) 
Recalling that M ts — E[a 2 s \Ft], a direct calculation shows that 



Ztr = 



M ts ds = A(T) + B{T)W t + C(T) (W 2 - t) 

+ D(T) (Wf - 3tW t ) + E(T) (W? - UW 2 



3r 



where the coefficients are 
A(T) 

B(T) 

C(T) 

D(T) 

E(T) 



a 2 (s) + sP 2 (s) + 



2 s{ s y< 



ds 



2(3{s)(a{s) + sS(s))ds 

(P 2 {s) + a(s)5{s) + 6 2 {s)s) ds 

6(s)/3(s)ds 



00 X 2 



ds 



Using again the fact that Vt = Z u , it follows that the state-price density is 

V t = A{t) + B(t)W t + C(t) (W 2 -t) 
D(t) {Wf-3tW t ) +E(t) (W t 4 -6tW t 2 +3t 2 ) . 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



(3.35) 
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Therefore bond prices in a one-variable third chaos model are given by the following 
ratio of forth degree polynomials in W t : 



_ A(T)+B(T)W t +C(T)(W?-t)+D(T)(W?-3tW t ) + E(T)(W?-6tW?+3t 2 ) 
" A(t)+B{t)Wt+C(t)(W?-t)+D(t)(W?-3tW t )+E(t)(W?-6tW?+3t 2 ) 



In particular, the initial term structure is given by 



A(T) It (" 2 (*) + a0*(a) + ^) 



ds 



p or = ^= — ~ , ; (3.37) 



ds 



3.4. Choices of the chaos coefficients. Inspired by the parametric forms 
used in descriptive models, we propose to model the function <j>\ itself in exponential- 
polynomial form. Observe that setting 

n ki 

Ms) =I>*(s)e- ClS , L 2 (s) (3.38) 

i=l j=0 

where Lq(s) = and c, > to guaranteed integrability, leads to a first chaos model in 
which forward rates are ratios of functions in exponential-polynomial form. We regard 
this rational exponential-polynomial family as a natural extension of the exponential- 
polynomial class considered in [T], with added flexibility and potentially better cal- 
ibration performance. In particular, this class allows us to reproduce all empirical 
shapes commonly observed for forward rate curves, such as increasing, decreasing and 
hump-shaped functions. 



Having made this choice for the function 4>i, equation (3.24) and its generaliza- 



tion (3.29) naturally lead us to consider functions a, /3,7 and S all belonging to the 
exponential-polynomial class as well. 

3.5. Calibration results. We describe the details of the term structure cal- 
ibration in Appendix |Aj We calibrate 14 different chaos models using two distinct 
data set from the UK bond market: observations of bonds of different maturities at 
every other business day from January 1998 to January 1999 (a volatile market often 
exhibiting an inverted yield curve) and weekly observations from December 2002 to 
December 2005 (a more moderate market). For comparison, we also calibrate two of 



the descriptive models specified in (3.1) 



We summarize our calibration results for the two data sets in Tables 13.11 and 13.21 
The first column in each table labels the models, starting with the descriptive Nelson- 



Siegel and Svensson models defined by the first two expression in (3.1), followed by 
the models (A.1)-(A.14) in Appendix [Aj The second column characterizes the type 



of model, whereas the third one gives the number of calibrated parameters. The 
remaining columns show the average values for the negative log-likelihood function, 
the Root-Mean-Squared Percentage Error and the Diebold-Mariano statistics with 
respect to the Svensson model, as described in Appendix |A"| 



Starting with the results in Table 3.1 for the volatile market, we first notice 



that, as expected, the Svensson model, which has 6 parameters, outperforms all three 
models that use a smaller number of parameters. On the other hand, out of five chaos 
models with 6 parameters, three are outperformed by the Svensson model, while the 
other two perform comparably. Finally, all chaos models with 7 parameters have 
a performance that is comparable to the Svensson model. We are led to conclude 
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Table 3.1 

Term structure calibration for 1998-1999 (Volatile Market) 





Model 


N 


-L 


RMSPE (%) 


DM 


Sv 


Svensson 


6 


160 


0.70 


~ 


NS 


Nelson-Siegel 


4 


2101 


2.67 


-4.45 


1 


1st chaos 


3 


4420 


4.44 


-11.46 


2 


1st chaos 


5 


250 


0.86 


-3.54 


3 


one-var 2nd chaos 


6 


162 


0.82 


-2.26 


4 


one-var 2nd cnaos 


7 


160 


0.69 


0.22 


5 


one-var 2nd chaos 


7 


145 


0.75 


-1.05 


6 


factorizablc 2nd chaos 


6 


335 


0.88 


-2.54 


7 


factorizablc 2nd chaos 


6 


245 


0.68 


0.27 


8 


factorizablc 2nd chaos 


6 


1245 


1.26 


-3.81 


9 


factorizablc 2nd chaos 


7 


179 


0.63 


1.38 


10 


factorizablc 2nd chaos 


7 


153 


0.72 


-1.07 


11 


one-var 3rd chaos 


6 


168 


0.72 


-1.24 


12 


one-var 3rd chaos 


7 


141 


0.76 


-1.16 


13 


one-var 3rd chaos 


7 


152 


0.72 


-1.19 


14 


one-var 3rd chaos 


7 


149 


0.76 


-1.43 



Table 3.2 

Term structure calibration for 2002-2005 (Moderate Market) 





Model 


N 


-L 


RMSPE (%) 


DM 


Sv 


Svensson 


6 


442 


0.76 




NS 


Nelson-Siegel 


4 


541 


0.97 


-1.76 


1 


1st chaos 


3 


8716 


3.96 


-3.50 


2 


1st chaos 


5 


438 


0.99 


-1.99 


3 


one-var 2nd chaos 


6 


388 


0.89 


-1.23 


4 


one-var 2nd chaos 


7 


388 


0.80 


-0.38 


5 


one-var 2nd chaos 


7 


329 


0.66 


1.26 


6 


factorizable 2nd chaos 


6 


437 


1.04 


-3.33 


7 


factorizable 2nd chaos 


6 


495 


0.84 


-0.68 


8 


factorizable 2nd chaos 


6 


421 


1.19 


-2.84 


9 


factorizable 2nd chaos 


7 


365 


0.82 


-0.78 


10 


factorizable 2nd chaos 


7 


323 


0.72 


0.36 


11 


one-var 3rd chaos 


6 


388 


0.87 


-1.06 


12 


one-var 3rd chaos 


7 


350 


0.78 


-0.11 


13 


one-var 3rd chaos 


7 


367 


0.68 


1.24 


14 


one-var 3rd chaos 


7 


325 


0.69 


0.60 



that chaos models do not offer a significant advantage with respect to the Svensson 
model under volatile market conditions, when yield curves can change from normal 
to inverted to humped shapes quiet rapidly. Moving on to Table |3.2| we observe 
smaller RMSPEs and higher DM statistics uniformly across all chaotic models, but 
nevertheless not high enough to reject the null hypotheses that they have the same 
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calibration performance as the Svensson model. 

More interestingly, both tables confirm our previous remarks in Section |3.3| that 
all chaos models behave similarly to a first chaos model from the point of view of 
calibrating the initial term structure. Indeed, we do not find significant difference in 
the calibration performances between chaos models of different orders. If anything, we 
notice that factorizable second chaos models perform slightly worse than their one- 
variable analogues, therefore justifying our assertion that one should avoid trying to 
calibrate the initial term structure using a complicated expression such as (3.23). 

We illustrate the results by plotting the RMSPE as a function of time in Figures 



3.1 and 3.2 For comparison, we restrict ourselves to the chaotic models with 6 



parameters in Tables [33] and [3T2[ as well as the Svensson and Nelson-Siegel models. 



One-var 2nd chaos 

Factrozable 2nd chaos (a) 
Factrozsble 2nd chaos (0) 



Factrozable 2nd chaos (c) 
One-var 3rd chaos 




Jan98 Mar98 Jun98 Sep98 Dec98 Feb99 

Calibration Date 



Fig. 3.1. Error for term structure calibration in 1998-1999. 

The basic overall conclusion of this section is that, as far as the daily calibration of 
initial term structure is concerned, we cannot state that there exists a significant dif- 
ference in the calibration performances between the Svensson and the chaotic models. 
However, as we mentioned before, chaotic models specify the stochastic evolution of 
interest rates, in addition to ensuring their positivity. In other words, the parameters 
calibrated from yield data may be applied directly to the volatility term structure. 
As we will see in the next section, this turns out to be greatly advantageous in the 
joint calibration of yields and option prices. 

4. Option calibration. In the previous section we showed that the chaotic and 
descriptive models have comparable performances when calibrated to yield curves. In 
this section, we take into account ATM European options, in particular caplets and 
swaptions. Since descriptive models for yield curves are silent regarding the future 
evolution of interest rates, we compare chaos models of different orders with other 
popular interest rate models in the literature, including short-rate models, models 
based on the potential approach and market models. But let us first analyze in 
greater detail the expressions for option prices in chaos models of different orders. 

4.1. Option prices in second chaos models. Since in a first chaos model the 
price at time of an interest rate derivative with payoff Ht reduces to the determin- 
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J One-var 2nd chaos 

6 - / Factrozable 2nd chaos (a) 

Factrozable 2nd chaos (b) 

— 1 — Factrozable 2nd chaos (c) 
— 3 — One-var 3rd chaos 



5 - _H -- Nelson-Spiegel 




Dec02 Jul03 Feb04 Oct04 May05 Dec05 

Calibration Date 



FlG. 3.2. Error for term structure calibration in 2002-2005. 



istic expression 



E\VtHt] Vt jj 

iJf) — — LIT 

° V V 



Io r - ds H 7 



(4.1) 



any nontrivial option pricing should start with at least a second chaos model. We 
begin with bond put options in a factorizable second chaos, for which it has already 
been observed in [IB] that (3.19) implies that 



KZ tt ~KZ tT =Vf\z), 

where z = Rt/y/Qt is a standard normal random variable and 

V { p 2) {z) := a + a x z + a 2 z 2 

is a second degree polynomial with coefficients 

a = -{A(T) - KA{t)) + (C(T) - KC(t))Q t 
ai = -(B(T)-KB(t))Q^ 2 (t) 
a 2 = -(C(T)-KC(t))Q t . 



(4.2) 



(4.3) 



(4.4) 



Using (2.14), it then follows that 

p(0,t,T,K) = 



1 



A{0)V2tt Jv^(z)>o 



~2 dz 



(4.5) 



As shown in |16) . once we calculate the roots of Vp (z), we have that (4.5) reduces 
to a simple expression given in terms of the standard normal cumulative distribution 
function. Having calculated the price of a put option on a bond, we can easily find 
the price of a caplet using expression (2.18). 
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Regarding the price of a payer swaption, it was also observed in [16] that 



Z tt - Z tTn - K^fr - T^)Z m = V™ p (z 



i=l 

where z is the standard normal random variable defined above and 

Vswp(z) := &o + hz + b 2 z 2 
is a second degree polynomial with coefficients 



(4.6) 



(4.7) 



bo = \A(t) - A(T n ) - K J2( T i ~ r <-iMPi)J 

- (c{t) - C(T n ) ~ K Y^ T i ~ T i-i) C i T ^ Qt 
bi = (B(t) - B(T n ) - K YS r ^ ~ Ti_ x )B{Ti)\ Q 1/2 t 
b 2 = [C(t) - C(T„) -K^T, - T i _ 1 )C(T i ) ) Q t . 



(4.8) 



It then follows from (2.24 1 that the price of a payer swaption in a factorizable second 
chaos model is 



Swp(0,T,A/\ K) 



1 



V^ p (z)e-^dz. 

,(2) 



A(0W2k Jvfl p {z)>o 

Once more, as shown in [16], having calculated the roots of V<f'(z), expression (4.9 1 



(4.9) 



reduces to an explicit formula in terms of the standard normal cumulative distribution 
and density functions. 

4.2. Option prices in third chaos models. We focus on one-variable third 
chaos models of the form (3.31). Recalling expression (3.331 for the process Z t T we 
have that 



KZ tt -Z tT = Vf\z) 
where z = Wt / \fi is a standard normal random variable and 

V^iz) := c + az + c 2 z 2 + c 3 z a + c 4 z 4 
is a forth degree polynomial with coefficients 



(4.10) 
(4.11) 



co = 


- (a(t) 




+ (p{T) - KC(tf) t - 3 (e(T) 


C\ = 


(b(t) 


- KB{t) 


t 1 ' 2 + 3 (3{T) - 5(t)) i 3/2 


C2 = 


(5(T) 


- KC(t f 


t + 6 (E(T) - KE(t)j t 2 


C3 = 


(p(T) 


- KD{t) 


) t 3/2 


c 4 = 


(e(t) 


- KE(t) 


t 2 



(4.12) 
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We then have the the price of a bond put option is a one-variable third chaos model 
is 

1 



p(0,t,T,K) 



A{0)V2tt Jv^(z)>q 



(4.13) 



which can then be used to obtain the price of a caplet according to (2.17). As before, 
once we calculate the roots of the polynomial Vp 4 \z) 7 both of these expressions reduce 
to explicit formulas in terms of the standard normal cumulative distribution and 
density functions (see [313]). 

Regarding the price of a payer swaption, we have that 



Z a - Z tTn - K^iTi - T^)Z m = Vfl p {z) 

where z — Wt/yt a standard normal random variable and 

^sl P ( z ) ■= d o + d x z + d 2 z 2 + d 3 z 3 + d 4 z 4 
is a forth degree polynomial with coefficients 

do = ^A(t) - A(T n ) - K YS r ^ ~ Ti-MiTi^j 

- \ C(t) - C(T n ) - K YS T * - Ti-^CiTM t 

+ 3 (S(t) - E{T n ) -KJ2( T * - Ti-i)E(T t )j t 2 
di = \S(t) - B(T n ) ~ K YS T * ~ Ti-iMT^ t 1 ' 2 

- 3 \D{t) - D{T n ) - K J2( T i - Ti-i)D(Ti)\ t 3 ' 2 
d 2 = (c{t) - C{T n ) - K YS r - - Ti-i)C(Ti)j t 

- 6 \E{t) - E(T n ) - K YS T * - Ti-iMT^ t 2 

d 3 = (5(t) - D(T n ) - K J2( T i ~ r<-i)5(Ti)^ t 3 

d 4 = [E(t) - E(T n ) -K^T, - Ti-jEpi)] t 2 



.3/2 



i=l 



(4.14) 



(4.15) 



(4.16) 



It then follows from (2.23) that the price of a payer swaption in a one-variable third 
chaos model is 



Swp{0,T,JV,K) = 



1 



-P { sl p {z)e-^dz. 



(4.17) 
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Once more, after we find the roots of ?1 4 ' (z), expression (4.17) reduces to an explicit 



formula in terms of the standard normal cumulative distribution and density functions 
(see [30]). 

4.3. Benchmark models. In this section we describe the models we use for 
comparing the calibration performance of chaotic models to option prices. 

We start with the simplest short-rate model for which the initial term structure 
can be freely specified in order to match observed data, namely the Hull- White [T7] 
dynamics 

dr t = - nr t )dt + r}dW t , (4.18) 
where we choose the function 6(t) in such a way that the initial term structure for 



forward rates has the Svensson form in (3.1 1, so that there are 8 parameters in total 



to be calibrated for this model. The prices for bond options, caplets and swaptions 
in a Hull- White model are given by well-known analytic expressions (see for example 
[1 pages 76-77]). 

Next we consider the rational lognormal model introduced in |llj (see related 
literature in [TJ], [20] and [53]). This is one of the earliest examples of a model in the 
potential approach and consists of setting 

a 2 s = gi {s)M s + g 2 (s), < s < t, (4.19) 

where gi , g 2 are nonnegative deterministic functions of time and M t is a strictly pos- 
itive continuous martingale such that M = 1. We then see that 

/>oo />oo 

Z tT = / M ts ds= / E[a 2 s \ F t ]ds = G 1 (T)M t + G 2 {T), (4.20) 

JT JT 

where 

POO POO 

Gi(t) := / gi(s)ds, and G 2 (t) := / g 2 (s)ds. (4.21) 



This implies that bond prices are represented in the following way: 

As shown in |llj . caplets and swaptions can be priced analytically in a rational log- 
normal model. Note however that, as also discussed in [TT] and [5], for each fixed t 
and T, bond prices in this model are bounded above and below by the determinis- 
tic ratios Gi(T)/Gi(i) and G 2 (T)/G 2 (t). As shown in [6], it follows that the short 
rate r t is also bounded above and below by the deterministic ratios G[ (t) / G\ (t) and 
G' 2 (t)/G 2 (t). As a consequence, rational lognormal models are unable to accurately 
price deep in the money and out of the money interest rate options. 

To implement the rational lognormal model, we follow [20j and parametrize the 
functions g\ and g 2 as 

gi (t) = -k 1 C ^(P ot ) k * and 32 (t) = _^[l- fcl (P 0t ) fc2 ], for t > 0, (4.22) 
for constants ki,k 2 6 R, where we choose the function P ot is then chosen so that the 



initial term structure of forward rates has the Svensson form in (3.1 ). Regarding the 



process M t , we follow 14J and represent is as an exponential martingale of the form 

M t = exp r)W t - ~T] 2 t , for t > 0, (4.23) 
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for some constant 77 € R. With this choices, there are 9 parameters in total to be 
calibrated in this model. 

The last benchmark that we consider are LIBOR market models. We first intro- 
duce the following notation: 

Fj(t) :=F m _ lTj , j = 1,2,..., 

where {To, T±, . . . , Tj_i, Tj, . . . , T n } is an increasing set of dates and F(t, S, T) is the 
forward LIBOR rate defined in ( |2.20 ). 

The lognormal forward market model (LFM), consists of specifying the dynamics: 

dF j (t)=* j (t)F j {t)dZ*(t), (4.24) 

where <jj(t),t > is a deterministic process, and Z J is a Brownian Motion under 
the forward measure Q- 7 , that is, the pricing measure associated with the numeraire 
P(',Tj) (see [I3j), with a correlation given by 

dZ l (t)dZ\t) = p i: jdt. (4.25) 

The motivation for introducing such models is that caplets can then be exactly priced 
using the Black formula introduced in (2.22). Following [2j page 224], we parametrize 
the volatilities as 

«7i(t) = h + (62 + 63(^-1 - t))^ 1 ^- 1 -*), (4.26) 

for 61,62, ^3, ci G R. Using the Svensson form in ( |3.1| for the initial term structure 
of forward rates, we see that there are 10 parameters to be calibrated to yield curves 
and caplet prices alone. 

To price a swaption maturing at t with tenor (T n — t) we use the Rebonato 
approximation (see [H page 283]) for the implied volatility 



v t,T n = 



\ 



T 



1 ^ ^(0)^(0)^(0)^(0) 



s(o,t,T n y 



<ji(t)aj{t)dt, 



(4.27) 



where the forward swap rates are assumed to be expressed as a linear combination of 
forward LIBOR rates as follows: 



S(t,T ,T n ) =J2^i(m(t) where w t (t) = {T ) Tl ' )PtT > 



(4.28) 



For the correlations appearing in (4.271, we use the Schoenmakers and Coffey [27] 
paramctrization: 



Pij = exp 



3 



n-1 



logPo 



Vi- 



i 2 + j 2 + ij — 3ni — 3nj + 3i + 3j + 2n 2 



(n-2)(rc-3) 



i 2 + j 2 + ij — ni — nj — 3i — 3j + 'in + 2 
(n-2)(n-3) 



(4.29) 



for i,j = 1,2, ...,n and constants 771,772,^00 G K satisfying the conditions 772 > 
and < T)x + 1]2 < — log poo. As we can see, if we use either only yields and caplets, 
we need to calibrate 10 parameters, whereas there are in total 13 parameters to be 
calibrated to yields and swaptions or yields, caplets, and swaptions together. 
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4.4. Calibration Results. We describe the details of the calibration to yields 
and option prices in Appendix [B] We calibrate 6 different chaos models using two 
distinct data set from the UK interest rate market: weekly observations of yields 
and implied volatilities for ATM caplets and ATM swaptions from September 2000 
to August 2001 and from May 2005 to May 2006. For comparison, we also calibrate 
the three benchmark models described in the Section l4~3l the Hull and White short 
rate model, the rational lognormal model and the lognormal forward LIBOR market 
model. 

We can gain a first impression in Figures 4.1 and 4.2 which show the calibrated 
implied volatility surfaces for ATM swaptions of different maturities and tenors at 
two specific calibration dates, one from each of the data sets above. We can see that 
neither the Hull and White nor the rational lognormal models have enough flexibility 
to reproduce the market quotes, with the LIBOR and chaos models providing much 
better fits. Figures |4.3| and |4.4| show similar results for the term structure of ATM 
caplets. In particular, we observe that both the Hull and White and the rational 
lognormal fail to capture the hump shaped term structure of caplet implied volatility, 
whereas the LIBOR and chaos models are able to do so. 

The quantitative calibration results are summarized results in Tables |4~l"1 to |4.6| 
The first column in each table labels the calibrated chaos models as listed in Appendix 
[5]and the benchmark models as listed in Sect ion l4~3l The second column characterizes 
the type of model, whereas the third one gives the number of calibrated parameters. 
The remaining columns in each table show the errors defined in Appendix [5] averaged 
across the different calibration dates. Specifically, the fourth column in each table 
shows the average value for the objective function minimized in the calibration. We 
plot this error across time for all calibration dates in Figures 4.5 to |4.10| The next 
columns show the average contribution to this total error coming from yields and the 
options used in each calibration. In the last column of Tables |4.1| |4.2| |4.4| and |4.5[ 
we also show the correspond to the pricing error obtained when we use the calibrated 
parameters to forecast the prices of options that were not used in the calibration. For 
example, in Table |4.1| we calibrate the model parameters to yields and caplets and 
use them to forecast the prices of swaptions. 

We can readily see that in general all chaos models that we consider fit the 
observed data better than the Hull and White and the rational lognormal models, 
as well as having smaller pricing errors for the options that were not used in each 
particular calibration. When compared with the LIBOR model, however, we see that 
only third chaos models have competitive performance. 

For a more detailed comparison, we consider relative frequencies for model selec- 
tion based on the Akaike Information Criterion described in Appendix [B] For each of 
our two data sets, we compared the relative performance of third chaos models with 
7 and 9 parameters and the benchmark LIBOR model. The results are presented in 
Tables 4.7 and 4.8 Whereas the relative performance of the third chaos model with 7 
parameters is somewhat ambiguous and depends on which type of options are consid- 
ered, we find that the third chaos model with 9 parameters consistently outperforms 
the lognormal forward LIBOR market model with 13 parameters. For example, in the 
first data set, it outperforms LIBOR in 36 out of the 53 calibration dates when we use 
yields and caplets, 53 out of 53 calibration dates when we use yields and swaptions, 
and 39 out of 53 calibration dates when we use yields, caplets and swaptions. For the 
second data set, the corresponding number of dates are 52, 44 and 39 respectively. 
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Table 4.1 

Yield and ATM caplet calibration for 2000-2001 



No. 


Model 


N 


TotalEl (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


5.1 


2.0 


4.6 


14.9 


2 


one-var 2nd chaos 


7 


3.3 


1.7 


2.7 


16.3 


3 


factorizablc 2nd 


6 


3.8 


2.1 


3.1 


26.5 


4 


one-var 3rd chaos 


6 


4.2 


2.0 


3.5 


15.5 


5 


one-var 3rd chaos 


7 


3.2 


1.3 


2.9 


15.7 


6 


one-var 3rd chaos 


9 


2.6 


1.1 


2.3 


17.0 


I 


Hull- White 


8 


8.7 


0.6 


8.7 


25.8 


II 


Rational-log 


9 


9.2 


0.6 


9.2 


13.9 


III 


LIBOR 


10 


3.0 


0.6 


3.0 





Table 4.2 

Yield and ATM swaption calibration for 2000 — 2001 



No. 


Model 


N 


TotalE2 (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


7.1 


1.8 


6.8 


14.5 


2 


one-var 2nd chaos 


7 


7.1 


2.0 


6.7 


14.6 


3 


factorizable 2nd 


6 


7.1 


2.1 


6.8 


14.3 


4 


one-var 3rd chaos 


6 


5.3 


2.9 


4.1 


10.2 


5 


one-var 3rd chaos 


7 


3.8 


1.5 


3.4 


8.6 


6 


one-var 3rd chaos 


9 


3.5 


1.5 


3.1 


9.1 


I 


Hull- White 


8 


10.2 


0.6 


10.2 


17.6 


II 


Rational-log 


9 


8.4 


0.6 


8.4 


15.3 


III 


LIBOR 


13 


5.0 


0.6 


5.0 


8.1 



Table 4.3 

Yield, ATM caplet and ATM swaption calibration for 2000 - 2001 



No. 


Model 


N 


TotalE2 (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


12.5 


2.2 


9.3 


7.9 


2 


one-var 2nd chaos 


7 


12.1 


2.4 


9.3 


7.3 


3 


factorizablc 2nd 


6 


12.1 


2.6 


8.4 


8.2 


4 


one-var 3rd chaos 


6 


8.2 


4.3 


4.4 


5.2 


5 


one-var 3rd chaos 


7 


7.1 


1.6 


4.4 


5.2 


6 


one-var 3rd chaos 


9 


5.9 


2.2 


4.1 


3.4 


I 


Hull- White 


8 


18.4 


0.6 


12.2 


13.7 


II 


Rational-log 


9 


14.6 


0.6 


10.0 


10.6 


III 


LIBOR 


13 


6.5 


0.6 


5.5 


3.1 
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Table 4.4 

Yield and ATM caplet calibration for 2005 - 2006 



No. 


Model 


N 


TotalE3 (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


6.3 


1.6 


6.1 


9.4 


2 


one-var 2nd chaos 


7 


3.4 


1.5 


3.0 


14.0 


3 


factorizable 2nd 


6 


4.3 


2.4 


3.4 


20.0 


4 


one-var 3rd chaos 


6 


4.9 


1.9 


4.4 


26.2 


5 


one-var 3rd chaos 


7 


3.6 


1.4 


3.2 


14.2 


6 


one-var 3rd chaos 


9 


3.3 


1.3 


3.0 


17.7 


I 


Hull-White 


8 


8.4 


0.4 


8.4 


16.3 


II 


Rational-log 


9 


9.4 


0.4 


9.4 


10.5 


III 


LIBOR 


10 


3.5 


0.4 


3.5 





Table 4.5 

Yield and ATM swaption calibration for 2005 — 2006 



No. 


Model 


N 


TotalEl (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


6.5 


3.2 


5.5 


32.2 


2 


one-var 2nd chaos 


7 


5.0 


1.5 


4.8 


11.9 


3 


factorizable 2nd 


6 


6.8 


2.3 


6.4 


13.7 


4 


one-var 3rd chaos 


6 


4.5 


2.2 


3.8 


21.2 


5 


one-var 3rd chaos 


7 


4.2 


1.6 


3.8 


13.4 


6 


one-var 3rd chaos 


9 


3.1 


1.3 


2.8 


14.2 


I 


Hull-White 


8 


9.5 


0.4 


9.5 


11.2 


II 


Rational-log 


9 


8.2 


0.4 


8.2 


10.9 


III 


LIBOR 


13 


3.5 


0.4 


3.5 


14.8 



Table 4.6 

Yield, ATM caplet and ATM swaption calibration for 2005 — 2006 



No. 


Model 


N 


TotalEl (%) 


YldE (%) 


CplE (%) 


SwpE (%) 


1 


one-var 2nd chaos 


6 


10.4 


2.5 


7.3 


6.9 


2 


one-var 2nd chaos 


7 


8.6 


1.5 


6.3 


5.6 


3 


factorizable 2nd 


6 


10.3 


1.9 


7.9 


6.2 


4 


one-var 3rd chaos 


6 


9.1 


3.3 


5.8 


6.1 


5 


one-var 3rd chaos 


7 


7.8 


1.8 


5.0 


5.5 


6 


one-var 3rd chaos 


9 


5.9 


1.4 


4.1 


4.0 


I 


Hull-White 


8 


14.0 


0.4 


10.1 


9.5 


II 


Rational-log 


9 


13.0 


0.4 


8.4 


9.9 


III 


LIBOR 


13 


6.2 


0.4 


4.8 


3.8 



5. Concluding remarks. We proposed and implemented a systematic way to 
calibrate chaotic models for interest rates to available term structure and option data. 
The calibration performance to initial term structures is comparable to that of tra- 
ditional descriptive models for forward rates with the same number of parameters, 
with the advantages of guaranteed positivity and consistency with a fully stochastic 
model for the time evolution of interest rates. When we include option data in the 
form of at-the- money caplets and swaptions, we see that chaos models perform sig- 
nificantly better than the benchmark Hull and White and rational lognormal models, 
and comparably to a LIBOR market model with a higher number of parameters. If 
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Table 4.7 

AIC model selection relative frequency (first dataset) 





Model 


Caplcts 


Swaptions 


Joint 


comparison 1 


one-var 3rd chaos (7 par) 
LIBOR (13 par) 


2/53 
51/53 


50/53 
3/53 


23/53 
30/53 


comparison 2 


one-var 3rd chaos (9 par) 
LIBOR (13 par) 


36/53 
17/53 


53/53 
0/53 


39/53 
14/53 


Table 4.8 

AIC model selection relative frequency (second dataset) 




Model 


Caplcts 


Swaptions 


Joint 


comparison 1 


one-var 3rd chaos (7 par) 
LIBOR 


14/53 
39/53 


23/53 
30/53 


7/53 
46/53 


comparison 2 


one-var 3rd chaos (9 par) 
LIBOR 


52/53 
1/53 


44/53 
9/53 


39/53 
14/53 



we take the number of parameters into account, a conservative information criterion 
shows that one of our chaos models consistently outperforms the benchmark LIBOR 
market model. 

The underlying reason for the superior performance of LIBOR market models 
when compared, for example, to Markovian short-rate models is the rich correlation 
structure they provide for LIBOR forward rates at different dates. Similarly, in chaos 
models the resulting short rate is in general not Markovian, and our calibration results 
show that an equally rich correlation structure can be achieved without having to 
model forward rates individually under each corresponding forward measure. 

The next step in our program is to consider options that are not at-the-money. 
The current industry standard for this is to use stochastic volatility models for either 
forward or swap rates. In a paper in preparation we show that chaos models naturally 
give rise to rates with stochastic volatility, and explore this fact to calibrate them to 
the smile observe in interest rate data. 

Expectedly, much work remains to be done, in particular on the interpretation of 
the financial interpretation of the chaos coefficients, where techniques such as principal 
components analysis might shed some additional light on the meaning of the most 
relevant calibrated parameters. Extensions to other classes of financial products, 
notably exchange rate derivatives, are also possible. We believe the analysis presented 
in this paper demonstrates the viability of chaos models as serious contenders for 
practical use in the financial industry and will stimulate further work in the area. 
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Appendix A. Term structure calibration procedure. 

We begin by listing our choices of chaos coefficients of different orders according 



to the general parametric form in (3.38). We start with first chaos models directly 



inspired by the Nelson-Siegel and Svensson parametric forms in (3.1), namely: 



<Ms) = (ki + Me~ ClS (A.l) 
Ms) = (h + b 2 s)e- ClS + b 3 se~ c * s . (A.2) 

Moving up one chaos order, we consider next one-variable second chaos models with 

a(s) = (h + b 2 s)e- ClS , p(s) = (b 3 + b 4 s)e- C2S (A.3) 
a(s) = he- ClS , (3(s) = (b 2 + b 3 s)e~ C2S + Me~ C3S (A.4) 
a{s) = fa + b 2 s)e~ ClS + b 3 se~ C2S , f3(s) = 6 4 se~ C3S (A.5) 



as well as factorizable second chaos models with 



a(s) 


= bie 


C1S , f3(s) 


= b 2 e 


C2S , 7(*) = 


(I + 63 


)e~ 


C3 S 


a(s) 


— b\e~ 


cis , m 


= (b 2 - 




7(5) = 


e~ 


C 3 5 


a(s) 


= (h- 


f b 2 s)e- ClS 


P{s 


= b 3 e- c * s , 


7(s) = 


e~ 


C35 


a(s) 


= he- 


C1S , f3(s) 


= (b 2 - 


\-b 3 s)e- c * s , 


7(s) = 


(1 


+ b 4 s)e- C3S 


a(s) 


= (h- 


¥b 2 s)e~ ClS 


P(a 


= b 3 e- c * s , 


70) = 


(1 


+ b 4 s)e~ C3S 



(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 



Finally, we consider the following one-variable third chaos models: 



a(s) 


— bie 


-Ci S 

1 


(3( S ) 


= &2e 


C2S , *(*) = 




(A.11) 


a(s) 


— bie' 


-Ci s 

) 


(3(s) 


= &2e~ 


C2S , 6(8) = 


(63 + &4s)e" C3S 


(A.12) 


a(s) 


= 6 x e- 


-Cl s 

) 


P(s) 


= (62- 




<J(s) = 6 4 e _C3S 


(A.13) 


a(s) 


= (61 


f b 2 s) 


e -cia 


/3( S 




<J(s) = 6 4 e- C3S 


(A.14) 



We now describe in detail the steps we take to calibrate these chaos models to 
observed yield curves, which we obtain from clean prices of treasury coupon strips in 
the UK bond market from the United Kingdom Debt Management Office (DMO) Q 
according to 

Vot ■= log Pot, (A.15) 

using an Actual/Actual day-count convention PJ. We consider the following two data 
sets: 

1. Bond prices at 146 dates (every other business day) from January 1998 to 
January 1999, with around 49 to 62 maturities for each date. 



http:/ /www.dmo.gov.uk/index.aspx?page=Gilts/Daily 
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2. Bond prices at 157 dates (every Friday) from December 2002 to December 
2005, with around 100 to 130 maturities for each date. 
Note that the first dataset contains a volatile market including the period of the Long- 
Term Capital Management (LTCM) crisis, whereas the second dataset corresponds 
to more moderate market conditions. 

We then apply the maximum likelihood estimation (MLE) method suggested 
by Cairns in [3] [7] to each of the models ( A.l )-( A.14[ ) and each starting date in 



the data sets above. This is done as follows: given a model with parameter vector 
= (b±, . . . , frjvj , ci, . . . , cjv 2 ) and a starting date to = with available maturities 
Tj, i = 1, ... ,71, we denote the theoretical prices by PoTi(0) and the corresponding 
observed bond prices by -PoT; ■ We then assume that 

log Pot, -AA(logPoT i (e),^ 2 (P T I (e),d l )), (A. 16) 

where di is the Macaulay duration for the bond and 



CT (P) d b (P) + 1 CT (P) Foo - <Jq (P)\ 

with the error parameters 



°o (p) = ooriri - ' ^ = 0.0005, tr 00 = 0.001, (A.18) 

adopted in [3] for the UK bond market data between January 1992 and November 
1996. As explained in [i-, cto(p) is based on the assumption that the published bond 
prices have rounding error of around 1/32 per 100 nominal price, whereas ad corre- 
sponds to the assumption that the difference between actual and expected yields have 
independent errors of the order of five basis points. Finally, (Too places a limit on the 
magnitude of price errors for long dated bonds. 
This leads to the log-likelihood function 

m = -l± [.o g[ ^(P.,(e),*)l + "^(e^'l ■ (*■»> 

We then use a global search procedure to find the global maximum for the log 



likelihood function (A.19), that is, to avoid finding a local maximum, we repeat the 
procedure using 1000 different random starting points and select the best maximiza- 
tion result. 

Having estimated the parameter vector 0, we denote by yoTi(©) the fitted yield 
for maturity Tj and by y 0T . the corresponding observed yield. We then define the 
fitting Root-Mean-Squared Percentage Error (RMSPE) as 



RMSPE = 



\ n t~i VoTi 



We then apply the Diebold-Mariano (DM) statistics [5] based on RMSPE to 
compare fitting performances as is done in [18] and |29j . Here for the computation we 
use the program DMARIANC0 in the statistics package STATA, with a lag order of 



2 http: / /ideas. repec.org/c/boc/bocode/s433001. html 
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thirteen in both of our two data sets. The null hypothesis, which is that two models 
have the same fitting errors, may be rejected at 5% level if the absolute value of the 
DM statistics is greater than 1.96. We compare the calibration performance of the 
Chaos Models with a descriptive model for forward rates in Svensson form: the higher 
the DM statistics, the more a chaotic model outperforms the Svensson model. 

Appendix B. Option calibration procedure. 

For option calibration, we consider one-variable second chaos models with 



a(a) = (h + Me~ ClS , = (6s + b 4 s)e' C2S (B.l) 

a(s) = (6 X + b 2 s)e' ClS + b 3 se~ C2S , (3(s) = 6 4 se" C3S , (B.2) 

factorizable second chaos models with 

a(s) = (&!+Me~ ClS , f3(s) = b 3 e~ C2S , j(s) = (B.3) 

and one-variable third chaos models with 

a(a) = &ie- ClS , p(s) = b 2 e- C2S , 5(a) = & 3 e~ C3S (B.4) 

a(a) = (bx + b 2 s)e- cis , f3(s) = b 3 e~ C2S : 5(a) = 6 4 e~ C3S (B.5) 



a(a) = (h + b 2 s)e^ s , f3(s) = (b 3 + 6 4 s)e" C2S ! , 5(a) = (b 5 + b 6S )e~ C3S , (B.6) 

as well as each of the three benchmark models described in Section [4~3l 

Regarding the data, zero-coupon yields are bootstrapped from the LIBOR, Future 
and Swap rates (see [5] for the detail of the bootstrapping technique) and interest rate 
option prices from ICAP (Garban Intercapital - London) and TTKL (Tullctt & Tokyo 
Liberty - London) via the Bloomberg Database^ We consider the following two data 
sets from the UK interest rate market: 

• Data between September 2000 and August 2001 at 53 dates (every Friday 
closing mid price) consisting of 

— zero-coupon yields with 17 maturities ranging from one month to 20 
years, 

— implied volatilities for ATM caplets with 37 maturities ranging from one 
to 10 years, 

— implied volatilities for ATM swaptions with 7 maturities ranging from 
one month to 5 years and 6 tenors ranging from one to 10 years. 

• Data between May 2005 and May 2006 at 53 dates (every Friday closing mid 
price) consisting of 

— zero-coupon yields with 22 maturities ranging from one month to 20 
years. 

— implied volatilities for ATM caplets with 77 maturities ranging from one 
to 20 years, 

— implied volatilities for ATM swaptions with 7 maturities ranging from 
one month to 5 years and 6 tenors ranging from one to 10 years. 

Note here that the option data corresponds to a part of the data in [29], where data 
was analyzed between August 1998 and January 2007. We obtain caplet implied 

3 Here we would like to acknowledge helpfulness of the Bloomberg help desk staff, who have aided 
greatly our understanding of the actual market data. We particularly wish to extend thanks David 
Culshaw, from ICAP, for his assistance. 
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volatilities by bootstrapping ATM cap implied volatilities observed in the market 
using the technique described in [IS], where the ATM caplets implied volatilities 
maturing at six months and nine months are obtained by constant extrapolation. 
The extrapolation is necessary to bootstrap the other ATM caplet implied volatilities, 
but when we calibrate the data, the extrapolated prices give us great errors. Hence, 
although we follow [12) and implement the extrapolation, we do not use those two 
short maturities for the calibration. Moreover, we observed some obvious outliers and 
corrected them accordingly. 

For each of the models and data sets above, we perform three distinct calibra- 
tions: first to yields and caplets, then to yields and swaptions, and finally to yields, 
caplets and swaptions. To define the objective function for each to these calibra- 
tions, denote the observed yields and prices of ATM caplets and ATM swaptions by 
y 0T . , Cpl(Tj, Katm) an d Swp(Tj, Katm)i their theoretical counterparties by yoT;: 
Cpl(Tj, Katm) an d Swp(Ti, Katm), and the corresponding mean square percentage 
errors by 



YieldE = 



1 m 

-E 



Vm\ - VoTj 

VoTi 



(B.7) 



CplE = 

and 



SwpE — 

For the calibration to yields and caplets, we then minimize the objective function 

TotalEl = \J (YieldE) 2 + (CplE) 2 . (B.10) 

Similarly, for the calibration to yields and swaptions, we minimize the objective func- 
tion 



TotalE2 = \J (YieldE) 2 + (SwpE) 2 . (B.ll) 

Finally, for the calibration to yields, caplets and swaptions we minimize the objective 
function 

TotalE3 = sj (YieldE) 2 + (CplE) 2 + (SwpE) 2 . (B.12) 

For each of these calibrations, we test the pricing performance of the calibrated mod- 
els. For example, after calibrating to yields and ATM swaptions, we use the model 
to price the ATM Caplets and compute the pricing error from market ATM Caplet 
prices. Whenever possible, for example in the Hull- White model, we calibrate initial 
yield curves and options separately by minimizing their respective square errors, since 
these involve distinct sets of parameters. Moreover, for the LFM model, we minimize 
the square error in implied volatilities rather than actual prices. 
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'C^T^Katm) - CpI(Tj, Katm) 



Cpl(T i: K 



(B. 



ATM J 



-, "3 

-E 



i=l 



Swp(T a , Katm) ~ Swp(r i; K A tm) 
Swp(Ti, Katm ) 



(B.9) 



To access the relative performance between the models, we use the Akaike In- 
formation Criterion (AIC) and the model selection relative frequency described in 
[3]. The AIC is formed with the maximized value of the likelihood function L for an 
estimated model and the number of parameters k in the following way: 

AIC = -2 log(L) + 2k. (B.13) 

For the least squares method under the assumed normality of residuals, this reduces 
to . 

AIC = 7i log ( — ) +2fc + nlog(27r) + ", (B.14) 
\ n J 2 

where RSS is the fitted residual sum of squares. Since nlog(27r) + 5 is a constant, we 
ignore this term and conclude that 

AIC = nlog(?^) +2k. (B.15) 

To define the model selection relative frequency, let us suppose we have two models 
and N calibration sets (for example the different calibration dates in our data). For a 
data set j G {1, • • • , N} we compute AIC denoted AlCj for one model and AIC^ 
for the other. Suppose the AIC of the first model is smaller than the AIC of the other 
model I times. Then the model selection relative frequencies (MSRF) for the first 
model and the second model are computed respectively by 

MSRF (1) = ~ and MSRF (2) = (B.16) 
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Fig. 4.2. Implied volatility for ATM swapti^jis on August 12th, 2005 (Blue: Market Quotes, 
Green: Theoretical Values). 
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Fig. 4.3. Implied volatility for ATM caplets on December 8th, 2000. 
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Fig. 4.4. Implied volatility for ATM caplets on February 10th, 2006. 
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Fig. 4.5. Total error for yields and ATM caplets in 2000-2001. 
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Fig. 4.6. Total error for yields and ATM swaptions in 2000-2001. 
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Fig. 4.7. Total error for yields, ATM caplets and ATM swaptions in 2000-2001. 
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Fig. 4.8. Total error for yields and ATM caplets in 2005-2006. 
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Fig. 4.9. Total error for yields and ATM swaptions in 2005-2006. 
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Fig. 4.10. Total error for yields, ATM caplets and ATM swaptions in 2005-2006. 



